Radiative collisional heating at the Doppler limit for laser-cooled magnesium atoms 
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We report Monte Carlo wave function simulation results on cold collisions between magnesium 
atoms in a strong red-detuned laser field. This is the normal situation e.g. in magneto-optical traps 
(MOT). The Doppler limit heating rate due to radiative collisions is calculated for 24 Mg atoms in a 
MOT based on the 1 So- 1 Pi atomic laser cooling transition. We find that radiative heating does not 
seem to affect the Doppler limit in this case. We also describe a channelling mechanism due to the 
missing Q branch in the excitation scheme, which could lead to a suppression of inelastic collisions, 
and find that this mechanism is not present in our simulation results due to the multistate character 
of the excitation process. 

PACS numbers: 42.50.Lc,42.50.Vk,2.70.Uu 



I. INTRODUCTION 

Laser cooling and trapping methods are an impor- 
tant ingredient in the recent achievements in the low- 
temperature physics of gaseous atoms and molecules 
In alkali atoms one can reach very low temperatures with 
Sisyphus and polarization gradient techniques, which sur- 
pass the Doppler cooling method in efficiency, and appear 
in magneto-optical traps without additional efforts 0- 
The drawback is that it becomes impossible to test the 
basic two-state Doppler cooling theory with these sys- 
tems. 

The appearance of sub-Doppler cooling is based on the 
hyperfine structure of the alkali atoms. The same hy- 
perfine structure is reflected in the complicated molecu- 
lar state structure of the quasimolecule formed by two 
colliding atoms. In dilute atomic gases binary collisions 
dominate the atomic interactions, and the quasimolecule 
states couple with the cooling and trapping light field. 
This leads to inelastic light-assisted collisions that cause 
either loss or heating of the trapped atoms Thus 
in alkali atoms the hyperfine structure leads to difficul- 
ties in testing the Doppler theory, as well as in modelling 
the collisional processes. 

The major isotopes of alkaline earth atoms, on the 
other hand, have no nuclear spin and thus the hyper- 
fine structure is missing. This makes the Doppler limit 
for temperatures the true laser cooling limit, and allows 
for testing the basic Doppler theory. Experimental stud- 
ies have been published recently for 88 Sr gg, and have 
been obtained for 24 Mg as well |7|. They indicate an 
intensity-dependent heating rate which prevents reach- 
ing the theoretical Doppler limit. 

Inelastic laser-assisted collisions are expected to de- 
pend strongly on intensity even when the atomic tran- 
sition is saturated (I s = 0.444 W/cm 2 for the ^o^Pi 
laser cooling transition of 24 Mg) 4]. In this region the 
problem is that most theoretical approaches are applica- 



ble only to the weak field situation, and can not handle 
the energy exchange when the energy change spectrum 
is continuous (as is the case for radiative collisional heat- 
ing). An important exception is the Monte Carlo wave 
function method (MCWF) dSI, 

which we apply here to 
the case of 24 Mg, to obtain the average energy increase 
per collision near the Doppler limit (at 1 mK for 24 Mg). 
Although simple in theory, the collisions between the al- 
kaline earth atoms at strong fields are complicated since 
the partial waves and rotational states of the ground and 
excited states (respectively) of the quasimolecule form a 
network of two infinite sets of coupled states ^(| . 

In this article we report the results of Monte Carlo 
simulations of collisional radiative heating, based on a 
truncated set of partial waves and rotational states. We 
derive a heating rate due to these collisions, and com- 
pare it to the photon scattering heating rate that sets 
the Doppler limit in the standard collisionlcss theory. 
We also outline a channelling mechanism which is pos- 
sible due to the specific structure and selection rules of 
the alkaline earth dimer quasimolecule, and which could 
lead to the suppression of inelastic processes, and test it 
with the simulations. Here we concentrate on the results 
provided by the simulations, and report the multitude of 
technical simulation details elsewhere. 



II. MAGNESIUM QUASIMOLECULE 

The lasers for cooling and trapping in magneto- 
optical traps (MOT) are normally detuned a few atomic 
linewidths (7 a t) below the atomic transition, which 
means that the attractive quasimolecule states have a 
finite internuclcar distance, the Condon point Rc, where 
they become resonant with the laser field. For 24 Mg and 
other alkaline earth atoms we have two such states, a 
strongly coupled state, and a weakly coupled 1 Ii g 
state [ill E3- 
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In the radiative heating mechanism, the colliding 
atoms on a partial wave I are excited at Rq to an attrac- 
tive quasimolecule state with rotational quantum number 
,/. The kinetic energy of the relative atomic motion in- 
creases as the atoms approach each other further, until 
they return to the ground state via spontaneous emission. 
The energy exchange due to such a process forms a con- 
tinuous distribution below the trap depth energy (above 
which the radiative heating becomes radiative escape). 

Since the 1 U g state is only weakly coupled to the light 
field, and provides mainly a channel for trap loss, we 
can concentrate in our modelling of heating on the 
state only (see Ref. for a more detailed discussion). 
The ground state partial waves and the excited rotational 
states are coupled by the selection rules I —> J — I ± 1, 
i.e., the Q branch for molecular excitation is missing. In 
the weak field limit one can consider each partial wave 
as paired independently to a rotational state, forming a 
set of independent two-state systems. At strong fields 
multiple couplings are possible, and one can move back 
and forth on the "angular momentum ladder" shown in 
Fig. 2] Note that for the bosonic 24 Mg atoms the ex- 
change symmetry allows only even values for I. 

J 12 3 4 
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/ 2 4 

FIG. 1: The coupling scheme for the ground state partial 
wave I and the excited state rotational state J. 

Channelling. One should note, as Fig. [21 illustrates, 
that each pair (Z, J) is coupled resonantly at a different 
value of Rc- Thus at modest field strengths one could 
consider a situation where the separation of the Condon 
points is valid while the excitation probability is close to 
unity. Due to the missing Q branch one can see that 
except for the s-wave, we always have a sequence / — ► 
J — I — 1 — ► Z 7 = J — 1 = 1 — 2 so that an incoming 
ground state partial wave I is channelled via two crossings 
into another ground state partial wave I" = I — 2, then 
reflected by the centrifugal barrier, returning the same 
route back to the original state Z. This process would 
reduce the collisions practically elastic (the s-wave still 
remains " unchannelled" , though) . 



III. COLLISIONAL HEATING RATE 

The energy increase from a single collision connects 
to a heating rate and the appropriate rate coefficient in 
the following way. The collisional heating rate kh(ti) de- 
scribes the kinetic energy change per unit volume and 
unit time. For identical particles it will be equal to 
T;Kffn 2 , where Kh is the rate coefficient for collisional 
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FIG. 2: Ground and excited state potentials for partial 
waves/rotational states up to I = 20 and J = 21. Open 
circles indicate the point of couplings forming a net of Con- 
don points. Here 8 = — 37 a t- The ground states have been 
shifted up by the energy of one photon. 

heating, and n is the atomic density. The factor of 1/2 
removes the doubling in collision counting for identical 
particles. The total collisional heating rate is obtained 
by integrating kh over the trap volume. 

For trap loss rate coefficient K\ oss one normally cal- 
culates the collision frequency times the probability for 
loss, i.e., wiossj where v is the relative velocity, and cri OS s 
is the cross-section for collisions, including the probabil- 
ity for a loss event to occur. This simple classical picture 
is connected to thermodynamics by assuming a thermal 
equilibrium distribution of velocities, f(v)dv, over which 
we take an average, and obtain K\ oss = {va\ oss ). The dy- 
namics of the two-body collision on the microscopic level 
then enters in calculating cti oss . 

To calculate the heating rate coefficient Kh is slightly 
more complicated, because in addition to the probabil- 
ity of an inelastic collision event to happen, we also need 
to estimate the amount of kinetic energy increase asso- 
ciated with it. Technically, we would have a continuous 
distribution of final energies corresponding to each ini- 
tial value of v, due to the randomness of the spontaneous 
emission events. For practical reasons we consider an 
averaged rate, i.e., we calculate the average change in 
kinetic energy per unit time. Since we can perform the 
averaging over final energy states before the averaging 
over initial states, wc define an inelastic (heating) cross 
section <th(v) (units energy x length 2 ) that gives the dif- 
ference of the average final relative kinetic energy and 
initial relative kinetic energy. 

In the partial wave approximation we can write the 
quantum mechanical cross-section for identical atoms in 
a three-dimensional trap as [l3| 

oo 

°W = W £ ( 2 * + 1 ) P ^)> (!) 

2=0(even) 
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where k is the wave number related to v and Pi{v) is the 
event probability. Thus, as a generalization, we write for 
heating 

oo 

aH ^ = ^ E (2l + l)^E sc (v,l), (2) 

/— O(cvcn) 

where AE sc (v, I) is the average single-collision energy in- 
crease related to the initial partial wave 

In principle one should take the thermal average over 
an isotropic Maxwell-Boltzmann distribution of relative 
velocities v. Due to the complexity of the Monte Carlo 
simulations we are, in practice, limited to calculating 
AE sc (v,l) for a rather narrow range of initial relative 
velocities v. As in Ref. we define a single-energy 

rate coefficient as (E — j^v 2 ) 

oo 

hQ E , ^ . 

with Q E = (2tt vE/h 2 ) 3 / 2 , where /i is the reduced two- 
particle mass. The simulations have shown this to be a 
reasonable approximation. 



IV. SIMULATION RESULTS 

We have performed a series of computationally tedious 
simulations at the detuning S = — 3j&t for various laser 
intensities, and the main results are given in Tabled Wc 
have used the semiclassical trajectory argument (turning 
point should be at R < Rc for any involved partial wave) 
to truncate the sum in Eq. © . Clearly a strong increase 
in the energy change with intensity continues above the 
saturation intensity I s . The semiclassical Landau-Zener 
theory also indicates that in this intensity region exci- 
tation probabilities should be equal to unity. Thus the 
continued increase can be attributed to the re-excitation 
of the decayed population, as suggested in Ref. 3; the 
delayed decay model predicts (roughly) AE SC oc yJTJT s . 



I (W/W) 


I/I, <V7at 


AE BC /E r 


0.036 


0.08 


-3.0 


16 ± 3 


0.88 


2.0 


-3.0 


145 ± 16 


2.2 


5.0 


-3.0 


249 ± 33 


3.6 


8.0 


-3.0 


307 ± 31 


5.3 


12.0 


-3.0 


394 ± 42 


8.0 


18.0 


-3.0 


427 ± 40 


14.2 


32.0 


-3.0 


547 ± 81 



TABLE I: Results from the multistate single-collision simula- 
tions. The total number of states has been 12, and the initial 
state has been 1 = 8. The energy increase AE SC is calculated 
as a time average of the kinetic energy in the region where it 
is flat, subtracted by the initial energy. The error refers to 
the statistical error of the average kinetic energy at the end 
of the Monte Carlo simulation. 



Additional simulations showed that pre- and post- 
collision steady state formation is fast and strong for 
I > 5/ s , and makes the energy exchange independent 
of the initial partial wave, i.e., AE sc (v,l) — > AE sc (v) in 
Eq. In other words, in a system where the energy 
states are strongly coupled and remain so even asymptot- 
ically, the definition of an initial state is obscured, and we 
can actually make meaningful studies because the spon- 
taneous emission always leads to the same steady state, 
independent of initial I, well before Rc is reached. For 
I < 5I S there is some initial state dependence, but on a 
scale that allows one to use the values given in Table [I] 
as a reasonable estimate of magnitude. 

As for the channelling effect, there was no sign of it 
in the simulations. The explanation is rather clear. The 
Condon points are close to each other both in position 
and energy, and the true picture of the excitation process 
is more like adiabatic following of the field-dressed (adi- 
abatic) quasimolecule states, shown in Fig. [21 However, 
this may not be the case for other values of detuning, 
or weaker fields than those used in our simulations. For 
this reason we want to emphasise the general possibility 
for channelling, and also because it is a good test for the 
validity of the reduction into independent two-state mod- 
els [l(| • We have confined our computationally demand- 
ing simulations to large intensities because then there 
are suitably many quantum jumps for having reasonable 
statistics for small ensembles of 64 or 128 members. 

Finally it should be noted that any truncation of the 
partial wave/rotational state manifolds must be done so 
that the total number of states is even. Otherwise one 
would introduce an artificial dark state as one of the field- 
dressed eigenstates, into which the strong and fast steady 
state formation would trap the system well before the 
actual collision begins. 




2it R/X 

FIG. 3: The field-dressed potentials for the magnesium quasi- 
molecule up to I = 12 and J = 13. Here 5 = — 37at and 
I/Is = 8.0. 
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V. COMPARISON OF RATES 

The heating rate prediction of the simulations must 
be compared to the photon scattering heating rate that 
enters the calculation of the Doppler limit. In this we 
follow Ref. and obtain the approximate rate 



dE 
~dt 



I/Is 



heat. 



1 + I/I s + (2V 7at ) 2 



(4) 



where E r — f?k 2 j(2\x) is the photon recoil energy, and 
I s = 7r7 a tft.c/(3A 3 ) is the saturation intensity 0. Adding 
the numbers for 24 Mg, we get the heating rate per particle 



dE 
~dt 



heat 



^ic60*. 
80 s 



(5) 



The above expression has been derived for a one- 
dimensional two-beam setup, and its extension to the 
full six-beam three-dimensional situation is an interest- 
ing and open question, and one of the strong motivations 
to study Doppler cooling with alkaline earth atoms. 

The collisional heating rate coefficient becomes, using 
Eq. ®, 



K H (E)=4n^%^(^ 
h \ E r 



1/2 



/3 7at \ 2/3 Air., 

\2\8\ 



(6) 



Here A = k~ x = A/(2tt), with A = 285.21 nm for 24 Mg. 
This result is obtained if we note that by truncating the 
partial wave series to the maximum classically allowed Z, 
denoted by Z max , we have E = h 2 l max (l max + l)/(2/xiJ 2 7 ). 
The Condon point Rc is determined by the 1 E+ state 
potential, U{R) = — 3h'f at /[2(k r R) 3 }, being equal to M 

at R c . Finally, I]'™o x (cvon) (2/ + 1) = (l max + l)(/max + 
2)/2 ~ / max (2ma X + l)/2 = [iR 2 c E / ' ti 2 . To compare with 
the single-particle rate JSJ) we write (for 6 = — 3 7a t) 



f) =\«»^ 

UL / coll Z 



fc B xl.3xl0" i4 — x 



K AE* C , 



Er 



■ xn. 

For n ~ 10 11 atoms/cm 3 and / ~ I s this gives about 0.1 
K/s. This result implies that radiative collisional heat- 
ing can not explain the problems in reaching the Doppler 
limit for 24 Mg, because the densities are usually in the 
range of 10 9 to 10 11 atoms/cm 3 . We have earlier per- 
formed two-state studies for Cs and found that there the 
heating rate present in Sisyphus cooling should match 
the collisional heating rate at n = 10 12 atoms/cm 3 |8j. 

As for making a comparison to the 88 Sr case, without 
performing simulations, we can make a rough estimate. 
In Eq. © we have written Kh(E) in a scaled form, so 
we can put in the corresponding values for 88 Sr as well 
(E/k B ^ T D = 0.77 mK, A = 460.73 nm, E r = 1 fjK, 
and taking again \S\ = 3 7a t)- The remaining problem 
is to estimate AE SC . In a toy model that assumes total 
excitation during a collision, and ignores re-excitation, 



we can write AE SC ~ U'(Rc)tv, where U'(Rc) is the 
slope of the excited state (ignoring J) potential at Rc, 
and r ~ 1/ 7at is the average survival time on the excited 
molecular state, and v is the initial collision velocity. This 
leads to an estimate of AE SC /E r = 23(E/E r ) 1 / 2 , which 
at least for 24 Mg gives 230, i.e., the correct magnitude 
near / — I s . With all this, we get for 88 Sr at the Doppler 
temperature the result Kh(E) ~ 2xl0~ 12 Kcm 3 /s, to be 
compared with E r -f at /80 ~ 2.5 K/s. These results meet 
only at densities of n ~ 10 12 atoms/cm 3 , which are a few 
magnitudes higher than those used in the experiment of 
Ref. This toy model also predicts a detuning depen- 
dence of Kh(E) of oc |6| 2 / 3 , i.e., heating decreases with 
detuning, which means that while the number of involved 
partial waves increases with diminishing \8\, the effect is 
more than countered by the reduction in the steepness of 
the excited state potential. 



VI. CONCLUSIONS 

This study is the first attempt to describe quantita- 
tively the radiative collisional heating by taking into ac- 
count the inherent multistate structure due to the par- 
tial waves and rotational states. The truncation of the 
partial wave manifold has required careful testing of var- 
ious possibilities. Although it appears that collisional 
radiative heating does not play a role in the thermody- 
namics of the trap at realistic densities for 24 Mg (and 
other alkaline earth atoms), our work does not rule it 
out completely. Although assumed to be irrelevant here, 
the weakly coupled 1 TI ff state may play a role, as well 
as the two quasimolecule states corresponding to the re- 
pulsive potentials. Experimentally the role of collisional 
processes can be distinguished from other processes by 
looking at the scaling with density; this has already been 
used to separate the cold collision contributions from the 
background collisions in trap loss. While we aim to im- 
prove our modelling and to use our heating results in 
estimating radiative trap loss, we also expect that more 
experimental data will become accessible in the near fu- 
ture. It should be pointed out that other effects may 
also lead to heating. The methods to study e.g. ra- 
diative trapping (reabsorption of scattered photons) are, 
however, very different from the ones employed here. 
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